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Abstract 

We use a Harrison transformation on solutions to the stationary axisymmetric Einstein equations 
to generate solutions of the Einstein-Maxwell equations. The case of hyperelliptic solutions to the 
Ernst equation is studied in detail. Analytic expressions for the metric and the multipole moments are 
obtained. As an example we consider the transformation of a family of counter-rotating dust disks. 

The resulting solutions can be interpreted as disks with currents and matter with a purely azimuthal 
pressure or as two streams of freely moving charged particles. We discuss interesting limiting cases 
as the extreme limit where the charge becomes identical to the mass, and the ultrarelativistic limit 
where the central redshift diverges. 

1 Introduction 

Electro-magnetic fields, especially magnetic fields are of astrophysical importance in the context of neu¬ 
tron stars, white dwarfs and galaxy formation. To get a complete relativistic understanding of such 
situations, one has to study the coupled Einstein-Maxwell equations. Since the stationary axisymmetric 
Einstein-Maxwell equations in vacuum in the form of Ernst [J| are completely integrable (see |2|), powerful 
solution generating techniques are at hand to obtain physically interesting solutions. But the equations 
in the matter region do not seem to be integrable, which makes it difficult to find global solutions which 
hold both in a three-dimensionally extended matter region and in vacuum. 

In the case of infinitesimally thin disks, the equations in the matter region reduce to ordinary differential 
equations the solutions of which determine boundary data for the vacuum equations. Such disks are 
discussed in astrophysics as models for Newtonian galaxies and the matter in accretion disks around 
black holes. Bicak, Lynden-Bell and Katz 0 studied static spacetimes, from which they removed a strip 
and glued the remainder together at the equatorial plane. The non-continuous normal derivatives at 
the equatorial plane lead to a (5-type energy-momentum tensor which can be interpreted as an infinitely 
extended disk made up of counter-rotating dust. This method was extended in to generate disk 
sources of the Kerr-metric. With the same techniques, disk sources for Kerr-Newman metrics static 
axisymmetric spacetimes with magnetic fields [B] and conformastationary metrics [7] were given. By 
construction all these disks have an infinite extension but finite mass. 

In the present paper we give solutions to the Einstein-Maxwell equations with disk sources of finite 
extension. The solutions are obtained by exploiting the SU(2, l)-invariance of the stationary vacuum 
Einstein-Alaxwell equations (see (S', to m- We use a so-called Harrison transformation m on disk 
solutions on hyperelliptic Riemann surfaces for the pure Einstein case, which are given in terms of 
Korotkin’s mi theta functional solutions, to obtain charged disks with currents. In this case analytic 
formulas for the complete metric can be obtained. As an example we discuss the transformation of a 
family of counter-rotating disks m which contains the Bardeen-Wagoner disk m (see for the exact 
solution) as a limiting case. Following |SJ the resulting disk can be interpreted either as a disk with a 
purely azimuthal pressure and currents or in a certain range of the physical parameters as two streams of 
freely moving charged particles, i.e. which move on electro-geodesics (solution to the geodesic equation 
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in the presence of a Lorentz force). We discuss the energy-momentum tensor at the disk, the multipole 
moments and interesting limiting cases. 

The paper is organized as follows. In section 2 we collect results on the stationary axisymmetric case 
in the absence of electromagnetic fields, solutions on hyperelliptic Riemann surfaces, and the family 
of disks discussed in 113 eg. The Einstein-Maxwell equations, their group structure and the action of 
Harrison transformations on the multipoles are discussed in section 3. In section 4 we consider a Harrison 
transformation of the counter-rotating disks m and discuss the metric and the energy-momentum tensor 
at the disk. The results are summarized in section 5 together with some concluding remarks. 


2 Stationary axisymmetric Einstein equations, theta functional 
solutions and counter-rotating dust disks 

In this section we will summarize results on the integrability of the Ernst equation, a class of solutions 
on hyperelliptic Riemann surfaces m and a member of this class describing counter-rotating dust disks 
which were discussed in m and CHI 


2.1 Ernst potential and metric 

We use the Weyl-Lewis-Papapetrou metric (see e.g. E?0] ) 

ds 2 = —e 2C/ (d t + ad(f>) 2 + e~ 2U (e 2fe (d p 2 + d( 2 ) + p 2 ^ 2 ) , (1) 


where p and £ are Weyl’s canonical coordinates and dt and d$ are the two commuting asymptotically 
timelike respectively spacelike Killing vectors. With £ = £ — ip, we define the complex Ernst potential 
£ = / + ib which is subject to the Ernst equation m 




(% - %) = 


£ + £ 


,£(£: 


tH 


( 2 ) 


where a bar denotes complex conjugation in C. If the Ernst potential is real, the corresponding spacetime 
is static and belongs to the so-called Weyl class. For a given Ernst potential the metric m follows from 
e 2U = f and 


a ? 


k 


2 P 


(g-g)e 

{£ + £) 2 ' 


(£-0 


(.£ + £) 2 ' 


(3) 

(4) 


Thus the complete metric can be obtained from a given Ernst potential via quadratures. 

Solutions with equatorial symmetry, i.e. a class of solutions where the metric functions have a reflection 
symmetry at the equatorial plane ( = 0 (this implies £(—() = £(()), are of special physical interest. 
In a Newtonian setting it can be proven that perfect fluids in thermodynamical equilibrium lead to 
equatorially symmetric situations, and the same is assumed to hold in a general relativistic context. A 
consequence of this condition is that Newman-Unti-Tambourini-parameters are ruled out. 

The stationary axisymmetric Einstein equations in vacuum were shown to be completely integrable in m 
and m- This implies that the Ernst equation has an infinite dimensional symmetry group, the Geroch 
group (see EH). The equation can be treated as the integrability condition for an overdetermined linear 
differential system for a 2 x 2-: 


'kcT ” 1 


d^-vp - 1 


itrix \P (see m) 



0 M 
0 


K - f / 0 A T 
K — £\ M 0 


(5) 
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here K is the so-called spectral parameter, and 


M = 


Jk_ 

5 + 5’ 


N = 


% 

£ + £' 


( 6 ) 


The spectral parameter is defined on the family of Riemann surfaces So given by the algebraic relation 
between and K, 

i4 = {K-£)(K-0. ( 7 ) 

Notice the special feature of the Ernst equation that the branch points £, £ depend on the spacetime 
coordinates. 

It is a consequence of the linear system © that a matrix T considered as a function on So and subject 
to the following conditions leads to a solution to the Ernst equation: 

I. The matrix U/ is at least twice differentiable with respect to £ and £. 'ff(P) is holomorphic and invertible 

at the branch points £ and £ such that the logarithmic derivative has a pole at £ and has 

a pole at £. 

II. All singularities of U/ on Eo are such that the logarithmic derivatives and \1/^R/ 1 are holomor¬ 

phic there. 

III. 'I' is subject to the reduction condition 


^{P G °) = a 3 V(P)a 2 , 


( 8 ) 


where 02 , 03 are Pauli matrices, and where <j 3 denotes the involution on Eo which interchanges the sheets. 

IV. The normalization and reality condition 


*1 '(P = oo + ) = 


£ —i 


(9) 


The function £ in © is a solution to the Ernst equation. 

We note that the choice of a 2 = ^ ^ ^ in condition III is a gauge condition, which does not fix the 

gauge completely, however. The remaining gauge freedom is due to matrices C(K) = ki(K) 1 + n 2 (K)a 2 , 
where the k, do not depend on £, £ and obey the asymptotic conditions «i(oo) = 1 and k 2 (oo) = 0. The 
matrices C act on T in the form H /C(K). Since it is a consequence of © that det ’3/ = F(K)e 2U 
where F(K) is independent of £, £ we can use this gauge freedom to choose F(K) = 1. The linear 
system © leads for the matrix y(P) = 'F _1 (oo _ )'I , (P) to one of the linear systems used in [2~T i. This 
parametrization, especially 


T 1 (oo )\I/(oo + ) 


1 


£+£ 


255 i(£ - £) 

%{£-£) 2 


( 10 ) 


reveals that the Ernst equation is invariant under 517(1,1)/50(1,1) transformations. For more details 
on the group aspect see m and the discussion in the next section. 

For given 'F(A"), one can directly determine the metric function a without having to integrate relation 
©. We denote by DpF(P) the coefficient of the linear term in the expansion of the function F in the 
local parameter near P. With © one finds for the matrix S = 1 (oo + )-D 00 +’P((X) + ) 




This implies with © 


($" 1 (c» + )D 00+ $) £ = 


a — a 0 = —25 


1 

( (55 _ ) c i(£ - £)t \ 

(11) 

2(5+ 5) 2 1 

[i(£%-£%) -(55) £ )■ 

to 

II 

to 

1 (oO + )P 00 + v I')i2, 

(12) 


where a 0 is a constant which is fixed by the condition that a vanishes on the regular part of the axis. 
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2.2 Solutions on hyperelliptic surfaces 

Using the integrability of the Ernst equation, Korotkin '5 gave solutions on hyperelliptic Riemann 
surfaces E g of genus g which are defined by the algebraic relation 

g 2 = (K - 0(K - 0 f[(K - Ei)(K - E t ). (13) 

i=l 


The points Ei = cti—iPi have to be constant with respect to the physical coordinates. We order the branch 
points with Im Ei < 0 in a way that $tEi < < ... < $lE g and assume for simplicity that the real parts 

of the Ei are all different. On this surface we introduce a canonical basis of cycles (a*, bi), i = 1,... ,g, 
the g differentials of the first kind d uji normalized by the condition j> d u>j = 2ni5ij , and the Abel map 

LOi(P) = Jp o d u>i. Furthermore, we define the Riemann matrix II with the elements 7= <f b d ojj, and 
the theta function 


0HW= ex p{^{ n ( iV + -^)>( 7V + -^)) + (( 2: + 7r * TO2 )’( Ar + f))} ( 14 ) 

NGZ3 1 ' < ' 


with half integer characteristic [to] = and to), to? = 0,1 ((N, z) = ^i z i)- We write O for the 
theta function with characteristic [0], and P = ( K , ±g(K)) or AT ± for a point PgE 9 (the + sheet will be 
defined by the condition that limi^,*, g(K)/K 9 = +1 and similarly for the —sheet). The normalized (all 
a-periods zero) differential of the third kind with poles at P\ and P 2 and residue +1 and —1 respectively 
will be denoted by d ujp 1 p 2 . 

Let lnG(AT) be a Holder-continuous function subject to the condition G(K) = G(K) and 


I=IJ In G(r )dui 00 + 00 - (r), m = J In Gdut 


(15) 


where T is the covering of the imaginary axis in the +-sheet of E g between —* and i. Then one obtains 
solutions to the Ernst equation in the form 


= 6[to](u(oo+) + u + u)(Q) j 
0[ra](u;(oo + ) — u + u;(£)) 


(16) 


Notice that these solutions depend only via the branch points of E g on the physical coordinates. In 
m it was shown that solutions of the above form on a Riemann surface of even genus g = 2s given by 
g 2 = (K — £)(A — £) UU(K 2 ~ 2 — Ef) with a function G subject to G(— AT) = G(K) lead to an 

equatorially symmetric Ernst potential. 

The Ernst potential <P|) follows with © and J = f r \n Gdujppcr from a matrix T of the form 


^ = e //2 / det(oo) / 

y det(AT) y 

with det(A') = 0[m](u + u;(P))0[m\(u — ui(P) + w(£)) + 0[m](zt — w(P))0[m](u + u>{P) + w(£)). 

To determine the metric function a via 1121) , we have to calculate the matrix S which leads with 11 7H to 


e[m](n+^(P)) J/ 2 
©[m] (u+o;(oo — ))_ 
e[m](M+i^(P)+u(g)) J! 2 
0[ra](u+a;(oo - )+a;(£)) 


■ e[m](«+o;(P g )) -J/2 

" 0[m] (u-\-cj(oo~ )} 

■ e[m](«+o;(P ,T )+^(g)) J/2 

0[ra] (u+cj(oo - )+o;(£)) 


(17) 


_ J_ 0(u + cu(oo )) 

12 2if °° + Q(u + uj(oo-)+u>(Oy 


and to determine the Harrison transformed function a in the next section, 


<?21 


££ 0(it + w(oo + )) 

2if °° + 0(u + w(oo+) + ui(£)) 


(18) 


(19) 
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Using a degenerated version of Fay’s trisecant identity mi (see m for the present case), we can write 
the above relations free of derivatives, 


p / 0[m](u)0[m](u + 2w(oo“) + u>(£)) 

2/ \L0[m](u + w(oo _ ) + w(£))0[to](u + u>(oo~)) 
p££ / 0[m](u)0[m](u + 2ol;(oo" , ') + w(£)) 

2/ \L0[m](zi + w(oo+) + w(£))0[m](u + w(oo+)) 


( 20 ) 


where 

r 0(^(OO~))0(^(OO~) +^(Q) , s 

0(O)9( W (O) ' ( j 

To construct the solution for the counter-rotating disks in [T^S], we used an algebraic approach which made 
it possible to establish algebraic relations between the metric functions at the disk. Let us recall that a 
divisor X on T, g is a formal symbol X = n\P\ + .. .+rikPk with Pi £ E and rii £ Z. The degree of a divisor 
is £T_i n-i■ The Riemann vector Kp is defined by the condition that 0(w(W) +Kr) = 0 if W is a divisor 
of degree g — 1 or less. We use here and in the following the notation uj{W) = fp dtu = Yli=i 
Note that the Riemann vector can be expressed through half-periods in the case of a hyperelliptic surface. 
We define the divisor X = Xu=i ^ as the solution of the Jacobi inversion problem (i = 1,... ,g) 


u>(X) — uj(D) = u, 


( 22 ) 


where the divisor D = J 2 i=i Ei (this corresponds to a choice of the characteristic [m] in ll 1 (ill 1. With the 
help of these divisors, we can write m in the form 


Inf 


r X 


ID 


t 9 dr 


-b [ i.cbt, 

27 ti J T p(r) 


(23) 


Additional information follows from the reality of u which leads to u>(X) — c o(D) = u(X) — u>{D). The 
reality condition for X implies via Abel’s theorem the existence of a meromorphic function R with poles 
in X + D and zeros in X + D (which is a rational function in the fundamental polygon), 


R(K) = const 


ULi( K - Ei)(K - £j) - Qo(K)p(K) 
UU(K - X Z )(K - Ei) 


(24) 


where Qq(K) = xq + X\K + ... + xK 9 1 is a polynomial in I\ with purely imaginary coefficients and 
x = ibe~ 2U . The coefficients a u are related to X via the relation 


9 9 

(1 - x 2 ) fj(A' - Xi){K - Xi) = fj(A- - Ei)(K - Si) - Q 2 0 (K)(K - £)(K - f). 

i= 1 2=1 


(25) 


We can use the existence of the rational function R to calculate certain integrals of the third kind as 


0(u + w(P))0(u + Lu(P a ) + w(£)) 
0(u + w(P <T ))0(u -(- ui(P) + w(()) 


rX+D 


= exp 


dupe 


p = 


Ix+D 


R(P) 

P(P CT ) 


(26) 


This makes it possible to give an algebraic expression for S12 and S21 ■ We can write S 12 in the form 


S 12 — 


1 


i{£ + £) 


A 


L ^ (L 


X+D 


fD 


I X+D 


diOpQ 


ID 


dl+pQ 


(27) 


with Q independent of P. The second integral can be reexpressed in terms of theta functions and be 
calculated with the help of so-called root functions: The quotient of two theta functions with the same 
argument but different characteristic is a root function which means that its square is a function on Y, g . 
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Let Pi, i = 1 ,... ,2g + 2, be the branch points of a hyperelliptic Riemann surface T, g of genus g and 
Aj = to(Pj) with w(Pi) = 0. Furthermore let {i\,... ,i g } and {ji, - • ■ ,j g } be two sets of numbers in 
{1, 2,... ,2 g + 2}. Then the following equality holds for an arbitrary point P £ S g , 


Q [Kr + ELi M P )) l (K-E il )...(K-E iB ) 

© [Kr + ELi J (w(P)) Cl V (K -E h )... (I< - E jg ) ’ 


where ci is a constant independent of K. 
This implies 

exp( / D 

Thus we get with 


S12 - ^ I 2 E( A "* ~ X i) + i ^ x 2 ( 2 ( ~ (£ + 0 

s. *=1 V \i=l y 


i(£ + £) 


and similarly 


<521 — ^ 


££ / 1 


-^)- t ~ e 2 f (+^)-(£+0 


i(£ + £) l 2 ■“ 
v 7 \ 2=1 


1 - z 2 l 2 




(28) 


(29) 


(30) 


(31) 


2.3 Counter-rotating dust disk 

We summarize results of m and m on disks which can be interpreted as two counter-rotating compo¬ 
nents of pressureless matter, so-called dust. The surface energy-momentum tensor S a/3 of these models, 
where a and (3 stand for the t, p and <f> components, is defined on the hypersurface £ = 0. The tensor S a P 
is related to the energy-momentum tensor T a @ which appears in the Einstein equations = 8ttT^ 1 ' 
(we use units in which the Newtonian gravitational constant, the dielectric constant and the velocity of 
light are equal to 1) via T a P = S a ^e k ~ u 5(Q. The tensor S a @ can be written in the form 


S a/3 = (J+U+U+ + a 


(32) 


where u± = (1,0, ±fi). We gave an explicit solution for disks with constant angular velocity fl and 
constant relative density 7 = ( a + — cr_)/(cr+ + cr_). This class of solutions is characterized by two real 
parameters A and <5 which are related to f 1 and 7 and the metric potential Uo = U (0, 0) at the center of 
the disk via 


A = 2n 2 e~ 2Uo 


(33) 


and 



(34) 


We put the radius po of the disk equal to 1 unless otherwise noted. Since the radius appears only in the 
combinations p/po, C/Po and ilpo in the physical quantities, it does not have an independent role. It is 
always possible to use it as a natural length scale unless it tends to 0 as in the case of the ultrarelativistic 
limit of the one component disk. The Ernst potential will be discussed in dependence of the parameters 
e = zr/(1 + zr) = 1 — e u ° and 7 , where zr is the redshift of photons emitted at the center of the disk 
and detected at infinity. 

The solution is given on a surface of genus 2 where the branch points of the Riemann surface are given 
by the relation Ei = —E 2 and E := E 2 = a + i/3 with a, (3 given by 



0 



(35) 
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The function G in m reads 


\J (t 2 - a) 2 + (3 2 + T 2 + 1 
” yj (r 2 — a) 2 + p 2 — (r 2 + 1)' 


We note that with a and 0 given, the Riemann surface is completely determined at a given point in the 
spacetime, i.e. for a given value of £. 

Regularity of the solutions in the exterior of the disk restricts the physical parameters to 0 < S < 6 S ( A) := 
2 ^1 + \J\ + 1/A 2 ^ and 0 < A < A c where A c ( 7 ) is the smallest value of A for which e = 1. The range of 
the physical parameters is restricted by the following limiting cases: 

Newtonian limit: e = 0 (A = 0), i.e. small velocities flpo and small redshifts in the disk. The function 
e 2U tends independently of 7 to 1 + A Un, where Un is the Maclaurin disk solution, and b is of order fl 3 . 
Ultrarelativistic limit: e = 1, i.e. diverging central redshift. For 7 7 ^ 1 it is reached for A c = 00 . The 
solution describes a disk of finite extension with diverging central redshift. For 7 = 1 , the limit is reached 
for A c = 4.629.... In this case the solution has a singular axis and is not asymptotically flat. This 
behavior can be interpreted as the limit of a vanishing disk radius. With this rescaling the solution in 
the exterior of the disk can be interpreted as the extreme Kerr solution (see m and references given 
therein). 

Static limit: 7 = 0 (6 = In this limit, the solution belongs to the Morgan and Morgan class m 

One component: 7 = 1 (<5 = 0), i.e. no counter-rotating matter in the disk. This is the disk of EH!, E3 
Analytic formulas for the complete metric in terms of theta functions are given in CHI. To evaluate the 
hyperelliptic integrals in the expressions for the metric we use the numerical methods and the cut-system 

of fT5] . 

At the disk the branch points £,£ lie on the contour T which implies that care has to be taken in 
the evaluation of the line integrals. The situation is however simplified by the equatorial symmetry 
of the solution which is reflected by the additional involution K —> —K of the Riemann surface £2 
for ( = 0. This makes it possible to perform the reduction K 2 —> r and to express the metric in 
terms of elliptic theta functions (see EH)- We denote with Y> w the elliptic Riemann surface defined by 
Vw = ( T + P 2 )(( r ~ a ) 2 + /® 2 )> an d let dw be the associated differential of the first kind with u w = 
A J _ 2 In G(y/r)dw(T). We cut the surface in a way that the a-cut is a closed contour in the upper sheet 
around the cut [—p 2 ,E] and that the 6 -cut starts at the cut [ 00 , .E]. The Abel map w is defined for 
P £ T, w as w(P) = f^ : dw. Then the real part of the Ernst potential at the disk can be written as 


e 


2(7 


1 

Y-S 




A 


+ 


+ s P 2 


(37) 


where 


y? + Sp 2 tf 2 (u w ) 
\J{p 2 + a) 2 + 0 2 VIM 


(38) 


It was shown that there exist algebraic relations between the real and imaginary parts of the Ernst 
potential, 





1 

A 




(39) 


and the function Z := (a — ao)e 2U 


Z 2 -p 2 + Se 4U 
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(40) 
















Moreover we have 


. _ z ( 1/X 2 + S l\ 

1X0 Sf [^i/x^ + Sp 2 a )■ ( } 

At the disk, the normal derivatives of the metric functions are discontinuous which leads via the Einstein 
equations to a 5-type energy-momentum tensor. To determine this tensor, it seems best to use Israel’s 
invariant junction conditions for matching spacetimes across non-null hypersurfaces l.'iOi . As in [191 we 
get 


si := -47rsoo 

s 2 := -47 t(s 0 3 - as 00 ) 

s 3 := -47r(s 33 - 2as 03 + a 2 s 0 o) 

where we define the tensor s^ via T^ = s 11 ” e 2 ( k ~ u ' > 8{(,). 
in terms of p-derivatives in the disk via 


= (k c ^2U ( )e 2U , 

= -\^ 2a , 

= -k c p 2 e~™, (42) 

The right-hand sides in can be expressed 


( 2 U\ Z 2 + p 2 + Se 4U 

(e )c = -—- b n 


k = - 


2 Zp 
P 2 


Z 2 + p 2 + 5e 4U , 


„2 U 


( 0 „ +V- 


2 Zp x Z 

This makes it directly possible to determine all quantities in the disk in terms of elliptic functions. 


(43) 


3 Einstein-Maxwell equations and Harrison transformations 

3.1 Field equations and group structure 

To treat the Einstein-Maxwell equations in vacuum in the stationary axisymmetric case we use again the 
metric CD We gauge the electromagnetic 4-potential in a way that it has only non-vanishing A t and 
components. It is subject to the Maxwell equations 

F tu p v = 0, *F^ V = ^e liva0 F a ^ = 0, F„ v = A^ v - A Vt „. (44) 

The Einstein equations lead to 

R pv = F^F? - (45) 

The vacuum Einstein-Maxwell equations in the stationary axisymmetric case are equivalent to the com¬ 
plex Ernst equations jQQJ 

fA£ = (V£ + 2$V<t>)V£, 

/AT = (V£ + 2<f>VT)VT, (46) 

where A and V are the standard differential operators in cylindrical coordinates, where the potentials 
£ and T are independent of </>, and where = / — TT and / = e 2U . The remaining metric functions 
follow via 

k = 

«C = 


8 


+ 2$$c)(f£ + 2T$e) - Me) 

J2 (i(^)£ + 


(47) 

(48) 







( 49 ) 


from £ and $. The electromagnetic potentials are related to $ via A t = and 


(^ = -f(A^-aA u ). 

For <f> = 0 one obtains the pure Einstein case of the previous section. 

The system is again integrable and can be treated as the integrability condition of a linear differential 
system as in the pure Einstein case (see EH). Even without axial symmetry the stationary Einstein- 
Maxwell equations are equivalent to a SU( 2,1)/£[Z7(1, 1 ) x [/(l)] nonlinear sigma model. This group can 
be parametrized by trigonal 3x3 matrices V , 

v7 


0 


0 


V = 


The matrix V satisfies 


iV 2 $ 1 0 

(b + m 2 )/Vf (V2i>)/VT W7 


V'rjV = rj, r) = 



(50) 


(51) 


i.e. it is unitary with respect to the metric r/ of SU(2 ,1). The action of Q £ SU(2,l) on V is V —> 
h(V,G)VQ~ 1 where h{V,G) £ 5[£/(l, 1) x XJ{ 1)] restores the triangular gauge of V. To obtain a gauge 
invariant parametrization, one introduces \ SWSE with 3 = diag(l, — 1,1) on which the action of 
Q £ SU(2, 1) is given by \ ^(G~ l )^xG~ l ■ We have 

If- 2|$| 2 + ( b 2 + |$| 4 )// V2$(6 - *|4>| 2 + if)/f (b - i\$\ 2 )/f \ 

X=\ -V2<S>(b + i\$\ 2 -if)/f 1-2|4>| 2 // -V2$/f • (52) 

V (b + i|$| 2 )// V2*/f 1 // / 

The SU( 2, 1) symmetry can be used to generate solutions by the action of an element G■ We list the 

0 0 0 

infinitesimal transformations and their consequences: 


add physically irrelevant constants to and 


changes / —> b, 



/ i6 

0 

0 > 

s. o 

0 

0 \ 

° 

- 2*0 

0 

is an electromagnetic duality transformation, 0 

0 

0 

V 0 

0 

id ) 

' V 0 

1 

o 


scale transformation, /, 6 , 4> 
which changes / —> $. 


e e f, e e b, e e / 2 4>, and 



are gauge transformations which 


is an Ehlers transformation which 


is a 


is a Harrison transformation m 


3.2 Harrison transformations 

Here we are interested to generate charged solutions from non-charged ones via a Harrison transformation, 
i.e. to generate solutions to the Einstein-Maxwell equations containing one additional constant parameter 
which is related to the charge. We assume that the original solutions are equatorially symmetric and 
asymptotically flat, i.e. / —> 1, <f> —> 0 and b —> 0 for \z\ —> oo. To ensure that the transformed solutions 
have the same asymptotic behavior, one has to use a scale transformation together with a transformation 
which changes 4> and b by some constant. By exponentiating the matrices of the SU (2,1) transformations, 
we consider a transformation of the form 


G = 


( 1 

id i 

-* 6 >i 0 i /2 \ / 


0 

0 \ 

0 

1 

- 0 i I diag(e _ ® 2 , 1 , e 6 * 2 ) 1 

-03 

1 

0 

V o 

0 

1 J \ 

— 04 + 0305/2 

-05 

1 / 


(53) 
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Since the asymptotic conditions imply that % and x! are the unit matrix at infinity, the matrix Q must 
satisfy the condition EQ'S = Q~ x . This leads with 15:111 to 


e * 2 = 


i-e^/ 2 1 


= iOi, 0 4 = 0, 


75 — 


The matrix Q thus takes the form 


1 _ id i_ —i0\0i/2 

G = - - TT To I ~ i - 1 1 + 6»i6»i/2 

- 0i0i/2 \ i6i Q i/2 __ 0i 1 


If we transform an Ernst potential in the case <f> = 0, we end up with 

, = 0i 0i0~i(/ 2 + & 2 )/2-(l + 0 1 0~ 1 /2)/+l-i6(l-0 1 0~ 1 /2) 

V2 (1 - 0i0"i//2) 2 + (0i0'i/2) 2 6 2 


(54) 


(55) 


(56) 


We are interested in transformations which preserve the equatorial symmetry, i.e. /(—£) = /(C), b(— C) = 
—6(C) and < f>(—C) = $(/). This implies for 15611 that 0i must be real which rules out magnetic monopoles. 
We put q = 0i/x/2 and sum up the results for the transformed potentials: 


= (1 - g 2 ) 2 / 

(1 — q 2 f) 2 + q 4 b 2 
(1 ~g 4 )0 
(1 — q 2 f) 2 + q A b 2 

. = _ (1 - /)(! - q 2 f) + q 2 b 2 + ib( 1 - q 2 ) 

q (1 — q 2 f) 2 + q A b 2 


(57) 

(58) 

(59) 


The real parameter q has to be in the region 0 < |g| < 1, for q > 1 the transformed spacetime would have 
a negative mass. The value q = 0 corresponds to the untransformed solution. The above formulas imply 
that the functions /', b' and 4>' are analytic where the original functions are analytic. 

The metric function k is invariant under the action of SU(2,1) transformations. To determine the 
transformed metric function a, we consider the matrix S in m- If we go over from 2x2 matrices to 
3x3 matrices according to the rule 


An A 12 

A 21 A 22 


An 0 A 12 

0 10 
A 21 0 A 22 


the matrix Q acts on S as on y. Thus we get with iTTTOl 


- (1 - fy 


(S 12 + 2,iq 2 — (? 4 52i), 


which is in accordance with BHll . This implies with GH> for the function a' 


(1 — 0 2 ) 2 


(5l2 — ? 4 <52l)- 


(60) 


(61) 


(62) 


To determine a' 0 , one has to consider £12 and £21 on the axis. In the limit p —> 0, there is a non¬ 
trivial contribution from the quotient of theta functions in (|20|) which diverges as 1/p. Repeating the 
considerations of m in the calculation of a 0 , one finds that the axis potentials can be expressed in terms 
of theta functions on the surface E of genus g — 1 given by p? = n?=i (k — Ei)(K — Ei). Denoting 
quantities on E with a tilde, one has with m on the axis 


S 21 0[m](w + 2<D(oo )) 
5i2 0[m](u + 2 cD(oo+)) 


(63) 
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In the equatorially symmetric case, this quotient is identical to one since 2d;(oo + ) is a half period on E 
(see EZD- Thus we have in this case 


l + g 4 
0 ~ 00 (1 — g 2 ) 2 ' 


(64) 


A well-known example is the Harrison transformation of the Schwarzschild solution which leads to the 
Reissner-Nordstrom solution. There we have with r± = -^/(£ ± m) 2 + p 2 the relations / = (r+ + r_ — 
2 to) / (r_|_ + r_ + 2m) ,6 = 0 which implies with 15011 


, = (r + +r_f-4m 2 = 2Q 

(r + + r_ + 2m') 2 ’ r + + r_ + 2m' ’ 


(65) 


and 6' = 0. We thus get a static solution with mass m' and charge Q which are subject to the relation 
m' 2 — Q 2 = m 2 . 


3.3 Asymptotic behavior 


We assume that the asymptotic behavior of the original solution, which is read off on the axis, is of the 
form / = 1 — 2M/|C|, 6 = —2 J/C, 2 and $ = Q/\C,\ — iJm/C, 2 plus terms of lower order in 1/|£| where M is 
the Arnowitt-Deser-Misner mass, J the angular momentum, Q the electric charge and Jm the magnetic 
moment. The same will hold for the Harrison transformed potentials. We find 


M' = M 


1 + g 2 


2g 


1 — g 2 1 — g- 


r Q, J' = J 


2q 


1 — g 2 1 — g- 


■ Jm, 


( 66 ) 


and 


Q' = Q 


l + g 2 _2g_ 

1 — g 2 1 — g- 


■M , J M — Jm 


l + g 2 _2g_ 

1 — g 2 1 — g : 


:J. 


(67) 


It is interesting to note that the quantities M 2 — Q 2 and J 2 — are invariants of the transformation. 
They are related to the Casimir operator of the 517(2, l)-group. If the original solution was uncharged, 
the extreme relation M' = ±Q' is only possible in the limit M = 0. 

A further invariant is the combination JmM — JQ which is of importance in relation to the gyromagnetic 
ratio 


2MJm 

9m = -jq~ 


( 68 ) 


Relation m implies that (Jm is equal to 2 if Q = Jm = 0 and q ^ 0. Thus all solutions which can 
be generated via a Harrison transformation from solutions with vanishing electromagnetic fields as the 
Kerr-Newman family from Kerr have a gyromagnetic ratio of 2. Due to the invariance of JmM — JQ 
under Harrison transformations, a gyromagnetic ratio of 2 is not changed under the transformation. 


4 Charged disks 

In this section we will use the Harrison transformation <1551 on the dust disk ns to generate a disk 
solution to the Einstein-Maxwell equation. We discuss the metric, several interesting limiting cases and 
the energy-momentum tensor. 

4.1 Metric functions 

Since the metric function e 2k is invariant under Harrison transformations, we refer the reader to llTij for 
a discussion. The metric function /' in C3 is proportional to /, which implies that the transformed 
solution vanishes exactly where the original solution has zeros. Since the set of zeros of / just defines the 
ergoregions, the transformed solution has the same ergoregions (if any) as the original solution. 
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For small q , the functions / and b are essentially unchanged since they are quadratic in q. The electro¬ 
magnetic potential $ is in this limit with lP|) of the form 

= — q(l — f + ib). (69) 

For larger g, \f'\ becomes smaller near the origin. Since its asymptotic values are not changed, the growth 
rate towards infinity increases which is reflected by the mass formula In the singular limit q —> 1, 

the function f' is zero for all finite values of z, but one at infinity. The behavior for b' is similar with 
the exception that b' is odd and zero at infinity. The function a also becomes singular in the limit q —*► 1 
which is reflected by the diverging factor 1/(1 — q 2 ) 2 and the constant ao (1SH which just implies that one 
can no longer choose a to be zero on the axis. In the metric function #03 = —a'/', the factors (1 — q 2 ) 2 
just cancel and the function is only marginally changed with increasing q. The typical behavior of g' 03 
for values of q with 0 < \q\ < 1 can be seen in Fig. Q] The metric function is an even function in £ which 


1.4.. . 

1 . 2 .. . 

1 

3.8 
a 

01 J.S 

3.4 

3.2 

0 
3 


Figure 1: Metric function g' os for e = 0 . 85 , 7 = 0.95 and q = 0 . 6 . 

vanishes on the axis and at infinity. It is analytic except on the disk where the normal derivatives have 
a jump. 

The electromagnetic potential tends to —1 in the limit q —> 1 for finite £, but is zero at infinity. The 
imaginary part is directly proportional to b' as can be seen from and <G2J). We show a typical 
situation for values of q with 0 < |g| < 1 in Fig. |2]for the real part and in Fig. 0for the imaginary 
part. The real part of $ is an even function in £ which vanishes at infinity and has discontinuous normal 
derivatives at the disk. The imaginary part of $ is an odd function in £ and has a jump at the disk. It 
vanishes at infinity. 

Since f has the same zeros as /, the transformed solution has a diverging central redshift if the un¬ 
transformed has, i.e. the ultrarelativistic limits coincide. In the case 7 / 1, one has a charged disk of 
finite extension with diverging central redshift. For 7 = 1, the solution in the exterior of the disk can 
be interpreted as an extreme Kerr-Newman metric which is obtained as a Harrison-transformed extreme 
Kerr metric. 

In the Newtonian limit A —> 0, one has / = 1 + ALbv and b = \ ?b in lowest order. This implies with <EJ 
to IHI for 1 — q 2 » A 

/' = 1 + b' = l±^\H = - 1 JL ^(XU N +i\%b), ( 70 ) 

1 — q z 1 — q z 1 — q z 
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Figure 2: Real part of the electromagnetic potential <f> for e = 0.85, 7 = 0.95 and q = 0.6. 


The transformed solution thus has the same Newtonian and post-Newtonian behavior as the original 
metric and in addition an electromagnetic field. The magnetic field is of order fl 3 as b'. 

Since the mass is of order A in the Newtonian limit, it is possible to have an extreme limit here with 
M' = Q' as in the Reissner-Nordstrom solution. If we put 1 — q 2 = reA with re > 0, we get in the limit 
A —♦ 0 for (BTfll and lt5§l) 


r = 


= 


Un 


(71) 


(re — U at ) 2 re — Un ’ 

a static solution similar to the extreme Reissner-Nordstrom solution, but with a jump in the normal 
derivatives of the metric functions at the disk and non-vanishing f 1 at the origin. Since Un < 0 in the 
whole spacetime, the solution is regular in the exterior of the disk. Thus one gets a non-singular limit in 
the exterior of the disk for q —> 1 in this case. 

In the static limit one has b' = = 0, since both are proportional to 5, 


/' 


(1 — q 2 f) 2 5 


= -q 


1 -/ 

1 -q 2 f 


(72) 


The Harrison-transformed static solution is thus again static with vanishing magnetic field but non-zero 
electric field. 


4.2 Disk and energy-momentum tensor 

It remains to be checked whether and in which range of the parameters the energy-momentum tensor 
m is physically acceptable. The above discussion of the metric indicates the extreme behavior of the 
metric functions for q close to one. It is plausible that the matter in the disk which is in the present 
example the source of such an extreme metric will in general not be physically acceptable. There can 
be maximal q smaller than 1 for given A and <5 which limits the physical range of the parameters. The 
matter in the disk will be interpreted by following Ledvinka et al. 0. 

To discuss the energy-momentum tensor and the currents in the disk, it is helpful to use the algebraic 
relations 1591 ) to 14111 between the metric functions which exist at the disk, and which imply similar 
relations between the transformed potentials. With m and ED, we get for S21 

S21 = £ SSi2 + ix 0 f, (73) 
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Figure 3: Imaginary part of the electromagnetic potential $ for e 


0.85, 7 = 0.95 and q = 0.6. 


and thus with for the metric function a' 

(1 _ q if(a - a 0 )' = (a- a 0 ) ^ + q ^ - jj + a + t) ) - ^ 

where a' 0 is given by m . With this function we can determine the angular velocity uji with respect to 
infinity of the locally non-rotating observers, for whom the metric is diagonal, 


u>i = - 


g 03 
533 


c 4 U 


p 2 — a 2 e 4C/ 


(75) 


This quantity is a measure for the dragging of the inertial frames with respect to infinity due to the 
rotating matter in the disk. The dependence of uii on e and 7 has been discussed in m■ As a function of 
q it is monotonically decreasing as can be seen in Fig. 0J The reason for this behavior is that the function 
f tends to zero in the limit q —* 1 for finite p, £ whereas g' 03 changes shape but remains finite. Thus 
the overall behavior of 07 is dominated by f. The deformation of the function g' 03 via q has however the 
consequence that u>i has its maximum for large q no longer at the center at the disk but near the rim. 
The energy-momentum tensor at the Harrison-transformed disk can be calculated via (TT2l) . Expressing 
the right-hand sides with the help of eo> and via the original functions, we get 

4 = ^ N 2 ^ ~ C 1 ~ g 4 / 2 + g 4&2 )/c + 2 g 4 /^c) » 

4 = 2^ ((! - 5 4 / 2 + q 4 b 2 )b p + 2 q%ff p ) , 

, _ p 3 N 

83 2(1 - q *ypz bp ' (76) 

where N = (1 — q 2 f) 2 + q 4 b 2 . 

Similarly a discontinuous electromagnetic tensor F leads to a current density J a via 

Ff = , , ( V-det (g^)F°A = -4ttJ“. (77) 

V-det (g M „) \ v /,/3 
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Figure 4: Angular velocity of the locally non-rotating observers with respect to infinity for e = 0.85, 
7 = 0.95, and q = 0,..., 0.9. 


Contributions to J arise only via the normal derivatives at the disk. We define the current density j a 
in the disk as s a @ by the relation J a =: e~ 2 ^ k ~ u ^ j a S(<^). The electromagnetic energy-momentum tensor 
does not produce a 5-type contribution to the Einstein-Maxwell equations since F is bounded at the disk. 
With (1771) we get using equatorial symmetry (the derivatives are taken at ( = 0 + ) 

2njo = 2?r(j 3 - aj 0 ) = -j(Ss®) p . (78) 

This implies 

2tt/o = q ^ N 2 9 ^ ("(C 1 ~ g 2 /) 2 ~ g 4&2 )/c + 2 g 2fe ( 1 ~ <ff)k) 

27r(j 3 - a jo) = ^ _ P g 2 j N f i 2 ^ 1 “ ^S)f P + ((1 ~ 9 2 /) 2 - Q 4 b 2 )b p ) . (79) 

The continuity equation (the Bianchi identity) at the disk T/"' = F IJV ,/„ = 0 leads to the condition 

3oo,pS°° + 23o3,pS° 3 + 533,pS 33 = 2(Fioj° + F\ 3 j 3 ). (80) 

We remark that one can substitute one of the equations 14211 by 18011 in the same way as one replaces 
one of the field equations by the covariant conservation of the energy momentum tensor in the case of 
three-dimensional perfect fluids. This makes it possible to eliminate k £ from iTTTll and to treat the energy- 
momentum tensor at the disk purely on the level of the Ernst equation. It is straight forward to check 
the consistency of this approach with the help of 03 • Thus one can solve boundary value problems of 
dust disks without using k. 

We will interpret the matter in the disk as in jj|>] in two ways. One possible approach is to introduce 
observers for whom the tensor s a @ is diagonal (</>-isotropic observers or FIOs) , i.e. we write in the 
form 

s a0 = a*v a v 0 +p*w a w 0 , (81) 
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where v and w are the unit timelike respectively spacelikc vectors N v (l,0,u^) and N w (k, 0,1) with 
v a w a = 0. We get 


_ ff33SQ0 ~ 90QS33 + \J (g335Q0 ~ gOQ533) 2 + 4(go3Spo — ffopSps) (g03S33 ~ ff33SP3) , , 

^ 2(g 0 3S33 — ff33Sp3) 


and 


go3 + ^<pg33 

goo + ^<l>g03 


(83) 


We show Ucj, for several values of q in Fig. 0 It can be seen that decreases monotonically with q. For 
large q the maximum of the angular velocity is near or at the rim of the disk in this example. 



Figure 5: Angular velocity of the FIOs with respect to infinity for e = 0.85, 7 = 0.95, and q = 0,..., 0.9. 


If the energy conditions are satisfied which is the case here, the FIOs can interpret the matter in the disk 
as a fluid with a purely azimuthal pressure. Alternatively they can interpret it as being made up of two 
streams of counter-rotating pressureless matter, S a ' 3 = \a* (t/“C/+ + !/“[/(!) with (U±) = iV(l,0, ±fi c ). 
The velocity f l c p is here below 1 (the speed of light). The matter streams will however not move on 
electro-geodesics unless j a w a = 0, i.e. if there are no currents in the frame of the FIOs. However it can 
be seen in FigEl that there are in general currents in the frame of the FIOs. Thus an interpretation of 
the matter by the FIOs as freely moving charged particles is not possible. 

A different approach is to interpret the matter in the disk as two streams of particles moving on electro¬ 
geodesics in the asymptotically non-rotating frame. To this end we make the ansatz 


s a ? = a+U$Ui + a-UlUi, j a = a+U? + *~U“ 


(84) 


with (U±) = N±(1,0,oj±). The angular velocity follows from the electro-geodesic equation for each 
component, 

< a- ...\ (85) 


? a m(9oo,p + 2gp 3 ,pW± + 533, pv±) = af (A 0 , p + A 3 , p uj±). 


Such an interpretation is possible if the angular velocities are real. Moreover the velocities uj±p in the disk 
should be smaller than 1 to avoid superluminal velocities, and the energy densities have to be positive, 
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Figure 6: Currents in the frame of the FIOs for e = 0.85, 7 = 0.95, and q = 0,..., 0.9. 


which means we are not interested in tachyonic or 


otherwise exotic matter. Relation (EH) leads to 


-Ni = 


■ W-J 


ui+ — w_ 


-Ni = 


u+r 


LU— — UJ+ 


( 86 ) 


r+Nl = 


s 00 

a~N 2 = - _ + 

KJ rrn -L V — « 


UJ-\. — UJ- 


J ^m ±y — 


UJ- — UJ-\- 


and 


UJ- = 


5 33 - LJ+S 03 

s 03 — u;_(_s 00 


If we enter ESI with this, we obtain 


(87) 

( 88 ) 


co± 


-T 2 ± ygf - FjTa 

7-i 


with 


,'0 „03 _ ,3 „00 

Ti = Oon - 2v4l ^ _A—_ 

1 1 933,p Z ^ l3,P S 03 s 03 _ s 00 s 33 

„-0 „03 _ .-3 „00 „-3 „03 _ a 0 „33 

— .903 .P ~ A °>P g03g03 _ s 00 s 33 ~~ A3 ’P s 03g03 _ s 00 s 33 

^03 j0 ^33 

7b — SOO.p — 2A 0iP ^ 03 ^ 03 _ s 00 s 33 ' 


(89) 


(90) 


The densities then follow from EH> where the continuity equation guarantees that this system can 
be solved. 

Numerically one finds that the angular velocities are real, but in a wide range of the parameters there 
are negative energy densities and tachyonic behavior. Already in the uncharged case there are infinite 
velocities in strongly relativistic settings with negligible counter-rotation which are due to extrema of the 
metric function (733 in the disk. In this case the quantity T\ in (p|) is zero which leads to a diverging 
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w_. Increasing q only enhances this effect. The result is that an interpretation as non-tachyonic counter¬ 
rotating matter on electro-geodesics with positive energy densities is only possible if q, e and 7 are not too 
large. In other words large values of q are in this setting only possible in post-Newtonian or nearly static 
situations. We show plots of the angular velocities u>± in Fig. [J for e = 0.36 and 7 = .08 where values 
of q up to 0.75 are possible. The corresponding densities are given in Fig. El and the charge densities in 
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Figure 7: Angular velocities u± for e = 0.36, 7 = 0.08, and <7 = 0,0.2,..., 0.8. 

Fig. El The densities vanish always at the rim of the disk, the charge densities are identically zero for 
q = 0. For larger values of q, the angular velocity becomes bigger and bigger at the rim of the disk 
and finally diverges. The density cr“ is in this case negative in the vicinity of the rim of the disk. 


5 Conclusion 

In this article we discussed Harrison transformations of solutions to the stationary axisymmetric Einstein 
equations in vacuum to generate solutions to the corresponding Einstein-Maxwell equations. The case 
of hyperelliptic solutions to the Ernst equation was studied in detail. This led to a class of solutions to 
the Einstein-Maxwell system which contains the following degrees of freedom: a real-valued function G, 
a set of complex branch points E, and a parameter q related to the total charge of the spacetime. Since 
the complete metric is given explicitly, one could in principle try to solve physically interesting boundary 
value problems as dust disks with a magnetic field within this class. 

As a first step in this direction we have considered a Harrison transformation of a family of counter¬ 
rotating dust disk solutions where the matter rotates in the frame of the FIOs on geodesics. The matter in 
the disk could be interpreted as moving on electro-geodesics in a certain range of the physical parameters, 
but in general not in the frame of the FIOs. As expected one cannot hope that a Harrison transformation 
preserves certain features of the matter as electro-geodesic motion in the FIO frame. If one wants 
to generate solutions with such properties or disks made up of charged pressureless particles without 
additional currents, one will have to solve the corresponding boundary value problem. Whether this 
will be possible within the class of Harrison transformed hyperelliptic solutions of the Ernst equation 
is an open question. For more general cases one will have to use Korotkin’s theta functional solutions 
to the Einstein-Maxwell equations on three-sheeted surfaces, where however the powerful calculus of 
hyperelliptic surfaces can no longer be used. 
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Figure 8: Matter densities a/ and er“ (dashed) for e = 0.36, 7 = 0.08, and q = 0,0.4, 0.8. 
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